sounding-base
Library to represent an atmospheric sounding with pressure as the vertical coordinate.
The base crate is meant to be a common base for other crates to build on. These crates may be for
managing a data-store, displaying data, or saving and loading files.
The emphasis of this crate is data representation and a common type for systems using sounding
data to build on and use.
Examples
use optional::{Optioned, some};
use metfor::{HectoPascal, Celsius, Feet};
use sounding_base::{Sounding, StationInfo};
fn main() {
let pressure_profile: Vec<Optioned<HectoPascal>> =
vec![1000.0, 925.0, 850.0, 700.0, 500.0, 300.0, 250.0, 100.0]
.into_iter()
.map(HectoPascal)
.map(some)
.collect();
let temperature_profile: Vec<Optioned<Celsius>> =
vec![13.0, 7.0, 5.0, -4.5, -20.6, -44.0, -52.0, -56.5]
.into_iter()
.map(Celsius)
.map(some)
.collect();
let stn = StationInfo::new_with_values(None, (45.6789, -115.6789), Feet(992.0));
let vt = chrono::NaiveDate::from_ymd(2018,3,8).and_hms(12,0,0);
let snd = Sounding::new()
.with_station_info(stn)
.with_valid_time(vt)
.with_lead_time(24) .with_pressure_profile(pressure_profile)
.with_temperature_profile(temperature_profile)
.with_station_pressure(some(HectoPascal(1013.25)))
.with_sfc_temperature(some(Celsius(15.0)));
let mut iter = snd.top_down();
let mut data_row = iter.next().unwrap();
assert_eq!(data_row.pressure, some(HectoPascal(100.0)));
assert_eq!(data_row.temperature, some(Celsius(-56.5)));
data_row = iter.next().unwrap();
assert_eq!(data_row.pressure, some(HectoPascal(250.0)));
assert_eq!(data_row.temperature, some(Celsius(-52.0)));
data_row = iter.next().unwrap();
assert_eq!(data_row.pressure, some(HectoPascal(300.0)));
assert_eq!(data_row.temperature, some(Celsius(-44.0)));
data_row = iter.next().unwrap();
assert_eq!(data_row.pressure, some(HectoPascal(500.0)));
assert_eq!(data_row.temperature, some(Celsius(-20.6)));
data_row = iter.next().unwrap();
assert_eq!(data_row.pressure, some(HectoPascal(700.0)));
assert_eq!(data_row.temperature, some(Celsius(-4.5)));
data_row = iter.next().unwrap();
assert_eq!(data_row.pressure, some(HectoPascal(850.0)));
assert_eq!(data_row.temperature, some(Celsius(5.0)));
data_row = iter.next().unwrap();
assert_eq!(data_row.pressure, some(HectoPascal(925.0)));
assert_eq!(data_row.temperature, some(Celsius(7.0)));
data_row = iter.next().unwrap();
assert_eq!(data_row.pressure, some(HectoPascal(1000.0)));
assert_eq!(data_row.temperature, some(Celsius(13.0)));
data_row = iter.next().unwrap();
assert_eq!(data_row.pressure, some(HectoPascal(1013.25)));
assert_eq!(data_row.temperature, some(Celsius(15.0)));
assert_eq!(iter.next(), None);
}
You probably noticed a lot of optional::Optioned
s in the example. Basically, anything can be
missing, and missing values are common in upper air soundings. For example, at high altitude the
dew point or humidity are often missing (if not totally inaccurate).